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Abstract 

We propose and study a model of polymer chains in a bilayer. Each chain is confined in one of the 
layers and polymer bonds on first neighbor edges in different layers interact. The thermodynamic 
properties of the model are studied in the grand-canonical formalism and both layers are considered 
to be Cayley trees. In the core region of the trees, which we may call a bilayer Bethe lattice, we find 
a very rich phase diagram in the parameter space defined by the two activities of monomers and the 
Boltzmann factor associated to the bond-bond interaction. Beside critical and coexistence surfaces, 
there are tricritical, bicritical end critical endpoint lines, as well as higher order multicritical points. 
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I. INTRODUCTION 



The theoretical study of the thermodynamic behavior of polymers, both in a melt or 
in solution, has a long history, continuous and lattice models were proposed to study such 
systems [l| • In lattice models, the first approach is to consider the linear polymeric chains to 
be random walks, this model is sometimes called ideal chain approximation. A more realistic 
approach is to represent the linear polymeric chains by self- and mutually avoiding walks 
(SAW's), and therefore one of the simplest models is athermal, so that only the excluded 
volume interactions are taken into account. Regarding the phase transitions which happen 
in such models, the excluded volume interactions are essential, at space dimensions below 
the upper critical dimension, to produce the correct critical exponents 2j. At and above 
the upper critical dimension the ideal chain exponents are found. In a grand-canonical 
ensemble, at low activities of monomers, a non-polymerized phase is stable, but as the 
activity is increased, a transition to a polymerized phase happens. This transition is usually 
continuous. 

While a simple model with excluded volume interactions may explain the polymerization 
transition when the chains are placed in a good solvent, in bad solvents and at low tem- 
peratures a collapse transition may happen. Below the transition temperature (0 point), 
the polymerization transition is discontinuous, so that the collapse transition temperature 
corresponds to a tricritical point. This transition is frequently named coil-globule transition 
in the literature. A simple effective model which explains this phenomena includes attrac- 
tive interactions between monomers placed at first neighbor sites which are not consecutive 
along a chain, these interactions favor more compact polymer configurations, thus reducing 
the contact area of the polymer chain with the solvent. In the literature, these interacting 
walks are called self-attractive self- avoiding walks (SASAW's) j^J. When it was realized by 
De Gennes that the simple polymerization model could be mapped on the ferromagnetic n- 
vector model in the limit when the number of components of the order parameter n vanishes 
jsj], the then new ideas of renormalization group, in the form of and expansion of critical 
exponents in e = 4 — d as a function of n, where d is the spacial dimensionality of the system, 
could be directly applied to models for polymers. Later, the relation between polymers and 
magnetic models was extended to SASAWs j^]. 

The behavior of magnetic models on bilayers and multilayers has attracted attention in 
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the literature for quite some time. For example, coupling through non magnetic layers shows 
oscillatory behavior and may lead to giant magnetoresistance effects jlj]. The phase diagrams 
of magnetically disordered bilayers have been studied on a system of two adjacent coupled 
Bethe lattices, for competingintralayer interaction parameters, the phase diagram displays 
multiple reentrant behavior [6j. A variety of other magnetic systems has also been studied 
on similar lattices [3] • Although, as expected, Bethe lattice solutions overestimate the region 
of stability of ordered phases, they may furnish the general features of phase diagrams. 

Polymeric chains close to interfaces between two immiscible liquids have also been studied, 
and it was found in a continuum model that the chains are attracted to the interface jsj, 
the maximum of the distribution of monomers being located close to the interface, inside 
the better solvent. A similar model was also studied in (q], with a different parametrization 
and several mean values and probabilities related to the conformation of the polymer chain 
are evaluated. The behavior of a copolymer close to a bilayer has been considered also 
in the literature, motivated by the relation of this problem with the relevant biological 
problem of the localization of a protein in a lipid bilayer 10[. Two regimes are found in 
this study regarding the localization of the chain in the bilayer. In the first, the density 
of monomers is a bimodal function, with maxima centered on the two interfaces, while the 
second configuration displays a single maximum of the density in the center of the bilayer. 
In the bimodal regime, few monomers are present in the area between the interfaces and 
most of the chain is adsorbed on them. 

Here we address a simpler problem, with homopolymer chains only, in the limit where each 
chain is entirely adsorbed on one interface and thus no polymer bond is present linking chains 
adsorbed on different interfaces. We include in the model interactions between polymer 
bonds placed in different layers on corresponding edges. This interaction takes into account 
the changes in statistical weight of a polymer bond if it is in contact with the other solvent 
across the interface or with a bond of another polymeric chain in the other interface. In figure 
CD a particular configuration of a part of a bilayer built with two square lattices is shown. 
In what follows, we consider both interfaces to be Cayley trees of arbitrary coordination 
number q, with polymeric chains on them whose endpoints are placed on the surface sites 
of the trees. Since we will study the behavior of the model in the core of the bilayer tree, 
we may consider this solution of the model to be a Bethe approximation for the model on 
regular lattices with the same coordination number. 
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FIG. 1: (color online) Part of a bilayer square lattice with a polymer chain on each layer. The 
three pairs of interacting bonds are edges of the shadowed (yellow) rectangles. 

The paper is organized as follows: In section [Til we define the model more precisely on the 
bilayer Bethe lattice and obtain its solution in terms of fixed points of recursion relations. 
The fixed points which are stable in some region of the parameter space are presented in 
section IIIIl studying the stability of these fixed points we address the transitions between 
the phases which are associated to each physical fixed point. The phase diagram if the 
three-dimensional parameter space of the model is studied in some detail in section \IV\ and 
final discussions and conclusion may be found in section [V] 

II. DEFINITION OF THE MODEL AND SOLUTION ON THE BILAYER BETHE 
LATTICE 

The problem is defined in the grand-canonical formalism, so that the parameters of the 
model are the fugacities of a polymer bond in layer 1, X\ = exp[ / u 1 /(A;sT)], and in layer 
2, X2 = exp^/^sT 1 )], as well as a Boltzmann factor x — exp[e/(A;BT)] for each pair of 
polymer bonds of corresponding edges in different layers, the energy of interaction of these 
pairs of bonds is, therefore, equal to — e, and is attractive if e > 0. In figure [2] a particular 
configuration of the model is shown for a bilayer Bethe lattice with coordination number 
equal to three. 

The grand-canonical partition function of the model may be written as 




(1) 
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FIG. 2: (color online) Bilayer Bethe lattice with q = 3 and M = 4 generations. The polymer 
configuration has a statistical weight x\ 4 x% X 3 = ^i 1 x \ K "\ where k = xix%X- 

where the sum is over configurations of self and mutually avoiding walks placed on the trees, 
N\ is the number of polymeric bonds on the Cayley tree 1, N 2 is the number of polymeric 
bonds on the Cayley tree 2, and Nb is the number of pairs of polymeric bonds located on 
first neighbor edges on different trees. We refer to figure [2] for the statistical weight of a 
particular configuration of the system. 

The solution of systems defined on tree-like structures defining recursive relations for 



partial partition functions (ppf's) is well known and is particularl y u seful for polymeric 



models because it is simple to impose the excluded- volume condition 12|, |_13|] . The general 



ization to a bilayer Bethe lattice is straightforward, and was widely used to study magnetic 
bilayers jg, 7|. We thus proceed defining partial partition functions for rooted sub-trees. The 
fixed root configuration of these sub-trees is specified by the polymer bonds which arrive 
at the root sites of sub-trees 1 and 2 coming from outer generations of sites. We need four 
partial partition functions: go,o; Qi,o] 9o,i an d <7i,i where the first (second) sub-index denotes 
the numbers of polymers bonds incident at the root site of the sub-tree 1 (2). Considering 
the operation of attaching q — 1 sub-trees to a new pair of root sites, we may write recursion 
relations for the partial partition functions of a sub-tree with one additional generation of 
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sites, denoted by primes. The recursion relations are: 

£0,0 = 00,0* + \{<1 - - 2 ) (&i0i,o + a$9o,i) 0o,o 3 + 
(g - 1) (g - 2) (g - 3) x 2 x 2 , x 0i,i 0i,o 0o,i 0o,o 4 + 
i(g-l)(g-2)x?x2 x 2 ( 7l 2 1 ^ 3 + 

i (g - 1) (g - 2) (g - 3) (g - 4) x 2 x 2 g 2 lfi gl , g^ 5 , (2a) 



0i,o = l) x i 9i,o 9o, 2 + (?-!)(? -2) 01,1 00.1 0o,o 3 + 

- (g - 1) (g - 2) (g - 3) x x x\ g lfl g 2 l , (2b) 



0o,i = (? ~ *) ^2 0o,i 0oV + (? -!)(?- 2 ) x 2 x 2 x 0i,i 0i,o 0o, O 3 + 

- (g - 1) (g - 2) (g - 3) a; 2 x 2 # 2 # ,i 0o,q 4 , (2c) 



0i,i = (9-1) 0i,i 0o,o 2 + (?-!)(?- 2) xix 2 g 1>0 00,100/ • ( 2 d) 

The ppf 's grow exponentially with the iterations, so, as usual, it is convenient to define 
ratios of the ppf's, and the thermodynamic properties of the model may be expressed by 
the fixed point values of these ratios. The ratios we use are: 



Ri = x 1 ^£ ; R 2 = x 2 9 -^ ■ R 3 = x 1 x 2X 9 -^ (3) 

00,0 00,0 00,0 

The recursion relations for the ratios will be 



(q - l)R x + (g - l)(g - 2) R 2 R 3 + -(g - l)(g - 2)(g - 3) i^i? 2 



R'r, 



(q - 1)R 2 + (g - l)(g - 2)R 1 R 3 + -(g - l)(g - 2)(g - 3) R\R 2 



X2 

D 



(4b) 



R' 3 = [(q-l)R 3 + (q-l)(q-2)R 1 R 2 ]-, 



(4c) 



where k = X\ x 2 x, and 
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D = 1 + i(g - l)(g - 2) + ^ + + ( g - - 2)(g - 3) R 1 R 2 R 3 + 

i(g-l)(g-2)( g -3)(g-4)^ (5) 

Connecting q sub-trees to the central site of the Cayley tree, the grand-canonical partition 
function takes the form 



p(q-l)(q-2)(q-3)R 2 1 R 2 2 . (6) 

The Bethe lattice solution of the model corresponds to its behavior in the core of the Cayley 
tree, so that surface effects are eliminated. The corresponding free energy per site may be 
found using a prescription proposed by Gujrati [l^]. The result is 

4> b = —7: [q In (D) - (g-2)ln(F)] (7) 



III. PHYSICAL FIXED POINTS 

The phase diagram is obtained looking for the physical fixed points of the recurrence 
relations Eqs.fllJ). For physical fixed point we understand that there exists a region in the 
thermodynamic space (xi, x 2 , k) where the fixed point is stable and corresponds to a global 
minimum of the free energy Eq.flT]), thus representing a thermodynamic phase. A fixed point 
R* = {R\,R^,R^} is stable if the Jacobian matrix 

dm 



dR j 



(8) 

R* 



has all the eigenvalues less than one. We found five stable fixed points describing five 
different thermodynamic phases: 
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1. Nonpolymerized fixed point (NP) 



The NP fixed point is given by r[ np) = 0; B%"'> = 0; R 



} (NP) 



? (NP) 



0, which corresponds to 



a phase where the density of polymers vanished in both trees. The Jacobian matrix is 



/ (9-1) 



J 



(NP) 



\ 



(q - 1) x 2 



\ o o 

thus the stability limits for the iVP phase are given by the conditions 

(NP) rwp\ 



X 



(NP) 



q-1 



x 



q-1 



-(NP) 



q-1 



(9) 



(10) 



2. Fixed point with polymers on the tree 1 (P±) 

The fixed point R\ ^ 0; R 2 = 0; Rs = 0, corresponds to a phase where the density of 
polymers is nonzero on the Bethe lattice numbered as one, and zero on the other lattice. 
The Jacobian matrix is given by 



/ 2-(g-l)xi 
(q-1) xi 



J 



(Pi) 



V 



o 



o 



[2+(g-l)(g-3)xi]x 2 
(g-1) X! 



2(q-2) [(g-l)xi-l] K 
q-1 xi 



For this phase, the stability limits are given by 







2(q-2)[(q-l)x 1 -l] X2 
q-1 Xi 



xi 



(11) 



X 



(Pi) 



q-1 



K — Xi 



2x 2 + (q - l) Xl [{q - 3)x 2 - 1] 
(q - l)[2x 2 - (q - l)xix 2 - xt] 



(12) 



3. Fixed point with polymers on the tree 2 (P 2 ) 



The fixed point R\ = 0; R 2 ^ 0; R3 = corresponds to a phase where the density of 
polymers is nonzero in the Bethe lattice numbered as two, and vanishing on the other lattice. 
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By symmetry, the stability limits of the P2 fixed point may be obtained from the stability 
limits of the Pi fixed point Eq. f JT2|) interchanging x\ and xi- 



4- Polymerized fixed point (P) 

For this fixed point R\ 7^ 0; R2 7^ 0; R3 7^ 0, so that the density of polymers is nonzero 
on both trees. In this case the stability lines were calculated numerically. 



5. Bilayer fixed point (Pb ) 

This fixed point Ri = 0; R2 = 0; R% 7^ corresponds to a phase where all polymer bonds 
are placed in pairs on corresponding edges of both trees, so that the polymer of both layers 
are fully correlated. The Jacobian matrix of this fixed point is 



/ 



j(P B ) 



2(g-2)((g-l) K -l) X2 
q—l K 



\ 
and the stability limits are given by 



2(g-2)((g-l)f t -l) xi q \ 

q—l k 



£2. 

K 



2-(q-l) K 
(9-1) « / 



(13) 



(q - - Xi)k 



(q - 1)(1 + 2(q - 2)x x )k - (3q - 5)x x 



k {Pb) 



(9-1) 



(14) 



IV. PHASE DIAGRAM 



The next step in order to obtain the phase diagram is to study the overlap between 
stability regions of different physical fixed points. If two regions where a unique fixed point 
is stable are separated by a common stability line, on this line a continuous phase transition is 
expected. If two or more fixed points are stable in the same region of the parameter space, 
this signals coexistence of phases. In this case, the loci of the discontinuous transitions 
may be found requiring that the coexisting phases should have the same (minimum) free 
energy. The five fixed points are present for coordination numbers q > 3, but for q = 3 the 
equations are simpler and thus some of the calculations may be done analytically. Therefore, 
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we performed all calculations below for q = 3. Some changes in the phase diagrams may 
arise for q > 3. 

Since a rather large number of physical fixed points is found, the model has a rich phase 
diagram, with different multicritical lines. In figures I3H we show six (£1,22) cuts for k = 
0; 1/4; 1/3; 9/20; 1/2; 1. These values of k were chosen in order to illustrate the main 
features of the phase diagram. The second-order lines and multicritical points were obtained 
analytically, and the first-order lines were found numerically. 

1 
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FIG. 3: (Color online) Phase diagram for k = 0. The dashed (black) curve is a first order transition 
and the full (black) lines are continuous transitions between the NP — P\ and NP — P2 phases. 
The bicritical point BCP is represented by a (violet) dot. This and the following diagrams were 
all obtained for q = 3. 



The first diagram, for k = 0, is depicted in figure [3j It actually corresponds to infinite 
repulsive interactions between polymer bonds in corresponding edges of both lattices. Al- 
though it should be said that this limiting case is probably quite far away from real polymeric 
systems, the phase diagram exhibits interesting features which lead us to shortly discuss it 
here. If a polymerized phase exists in one of the lattices, its presence inhibits the polymer- 
ization transition on the other lattice, since it reduces the effective coordination number at 
this other lattice. This effect is apparent in the diagram, the transition to the phase P2, at 
constant x% > l/(q — 1) happens at £2 = £1 > l/(^ — 1) and is discontinuous, since when 



10 



the free energy of the P± phase becomes smaller than the one of the P2 phase, x± is larger 
than its critical value l/(q — 1). The coexistence line of the two polymerized phases meets 
the two critical polymerization lines at a point which is a bicritical point (BCP), located at 
X\ — x-i — l/(q— 1). It is worth mentioning that the features of the transition lines incident 
at the BCP are the ones expected for a mean-field approximation as the one we are doing 
here: the classical crossover exponent is = 1, so that the critical lines are linear functions 



close to the BCP 



15l | . In general, one finds > and therefore the critical lines meet the 



coexistence line tangentially, as may be found, for instance, from e = 4 — d expansions for 



the ra-vector model when n > up to order e 4 



16] . However, all these coefficients vanish for 



n — > 0, so that it is possible that even below the upper critical dimension linear incidence 
of the critical lines at the BCP is observed. 




0,2 0,4 0,6 0,8 

x i 



FIG. 4: (Color online) Phase diagram for k = 1/4. The full (black) lines are continuous transitions. 
The green dot is a tetracritical point. In this and the following diagrams, on the dot-dashed (black) 
curve the interactions between the polymer bonds vanishes (x = !)• 



When k has a small positive value, the P phase occupies a region between the P\ and 
P2 phases, separated from them by critical lines, as may be seen in figure HI With this 
intermediate phase present, the bicritical point becomes a tetracritical point, similar to 



what is found in anisotropic antiferromagnets 



15]. Again the remarks about the crossover 



exponent <fi and the incidence of the four critical lines at the tetracritical point apply, and 
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for a mean-field approximation linear behavior is expected, in agreement with the results we 
obtained. In this and the following diagrams, the dot-dashed curve X1X2 = k separates the 
regimes of repulsive (x < 1) and attractive (x > 1) interactions between the polymers; the 
tetracritical point is located on this curve. The region of attractive bond-bond interactions 
if situated below the curve. In this phase diagram, the two critical lines which limit the P 
phase meet at a right angle. Between k = 1/4 and n = 1/3, this angle is larger than ir/2, as 
may be seen in the phase diagram for k — 1/3 shown in figure [5], where the angle is equal 
to 7r. A point which may be worth discussing is that in this case a reentrant behavior is 
seen in the critical lines which limit the phase P. This behavior is actually a consequence 
of our choice of variables, since k = X1X2X = const, implies that the Boltzmann factor of 
the interactions x is n °t constant if one of the activities changes. In a similar diagram with 
constant x t ne critical curves are monotonous and no reentrance is seen. 
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FIG. 5: (Color online) Phase diagram for k = 1/3. All transitions are continuous. The point where 
the four critical lines meet is the endpoint of a line of tetracritical points for k < 1/3 and four other 
lines, two of them tricritical and two of critical endpoints, for k > 1/3, thus being a multicritical 
point of higher order. The dot-dashed curve corresponds to x = 1. 

As n exceeds 1/3, a qualitative change of the phase diagram happens. The tetracritical 
point splits into two pairs of tricritical and critical endpoints, as may be seen in figure |6j 
A first order transition line now separates phase P and NP, while part of the transition 
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FIG. 6: (Color online) Phase diagram for k = 9/20. Two pairs of tricritical points (orange) and 
critical endpoints (blue) appear. The dot-dashed curve corresponds to % = 1. 

between phase P and the phases Pi and P2 becomes discontinuous. As k grows, the tricritical 
points approach the NP — Pi and the NP — P 2 critical lines and the critical endpoints, so 
that finally, at k — 1/2, the tricritical points meet the critical surface which limit the NP 
phase, as may be seen in figure [7J The multicritical points where the tricritical and critical 
endpoint lines meet are labeled MCP1 in the diagram. The dashed line which connects the 
two multicritical points is a line of bicritical points, since at this line two critical surfaces 
(NP — PB and PB — P) and a coexistence surface (NP — P) meet. At this value of k, the 
critical lines between polymerized phases change from concave to convex, so that they are 
linear. The NP phase is still stable in a region of the phase diagram close to the origin, but 
its stability is quadratic, so that the region covered by it in this phase diagram is actually a 
critical surface separating it from the bilayer polymerized phase Pb, which becomes stable 
for k > 1/2. On the lines between the multicritical points (MCP1 and MCP2), four critical 
surfaces meet (NP — P 2 , P2 — P, P — Pb an d Pb — NP on the upper part of the diagram), 
so that on this line we have tetracritical transitions. This feature may be better appreciated 
in the diagram in the variables k and X2, for x\ = 1/3, which cuts this line of tetracritical 
points and is shown in figure El 

As soon as k exceeds 1/2 all transitions are continuous. The phase diagrams are similar 
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x l 

FIG. 7: (Color online) Phase diagram for k = 1/2. The dotted (violet) curve is a line of bicritical 
points, and is limited by two multicritical points (red dots). Between these multicritical points and 
the multicritical points represented by squares, the critical lines P2 — P and P — NP, as well as 
Pi — P and P — NP coalesce, so that on this line four critical surfaces meet, thus characterizing 
it as a tetracritical line (green). The dot-dashed curve corresponds to % = 1. 

to the one shown in figure for k = 1. All four polymerized phases are present, there 
are three critical lines separating the phase P from another polymerized phase. Since the 
stability limits of the phases Pi, P2, and Pb are known (Eqs. (1121) and (j!4p ). the phase 
diagrams may be obtained analytically. The critical lines which limit the phases Pi and P2 
start at xi = and 22 = respectively, thus k = 1/2 is a quite singular value regarding the 
behavior of these lines. As expected, the region of the parameter space where the Pb phase 
is stable grows with increasing values of k. 

A. Identical Polymers 

It is interesting to study the case where both polymers are identical, that is x\ = £2 = %■ 
Besides the physical interest, this case is simpler than the general one. Only the phases 
NP, P, and Pb are present in this subspace of the general parameter space discussed above. 
Since in this case there is no particular advantage in using the composed variable k, we use 



14 



* 2 0,5 



0, 



1 1 1 .' 1 1 1 / 1 1 1 

/ / 


/ / 


/ / 


' p_ / 


/ 2 


/ p 


/ / 




1 1 




i 




i \ 




i 




i 




i 




i 




i 




i NP 




/ 




/ 




/ 




'.I.I 


I.I. 


0,2 0,4 


0,6 0,8 



K 



FIG. 8: (Color online) Phase diagram in the variables (re, X2) for x\ = 1/3. Four critical lines meet 
at a tetracritical point, represented by a circle (green). 




FIG. 9: Phase diagram for k = 1. The four phases present in this diagram are separated by 
second order lines. The dot-dashed curve corresponds to x = 1. 

the Boltzmann factor % instead. The phase diagram for q — 3 if shown in figure [TQ1 The 
line x = l/(q — 1) is a second-order line between phases NP and P from \ = up to a 
tricritical point localized at 
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x iTCX) = _L_ . X (TC1) = i+1 (15) 

For larger values of x, the -/VP — P transition is discontinuous. This part of the phase 
diagram is qualitatively identical to what if found for SASAW's, the tricritical point (TCP1) 
corresponds to the theta point where the SASAW's collapse. The NP — P coexistence line 
ends at a critical endpoint, whose localization has been determined numerically, since it 
involves the locus of the coexistence line. The phases Pb and P coexist on a line which 
starts at the critical endpoint and ends at the a second tricritical point (TCP2). Although 
it is possible to obtain an analytical expression for this tricritical point as a function of q 
using symbolic algebra programs, the resulting expressions are too long to be displayed here. 
In the particular case q = 3, the tricritical point is given by (xtc2 = (3v / 5 — 5)/4,xtc2 = 
16(5 — 2\/5)/(3\/5 — 5) 2 ). For values of x > Xtc2, the Pb — P line is a second order line 
given by 

X = i (l + (q - 2)x ± [(q~l)x(2+(q-2)x)-2}^ . (16) 

An interesting reentrant behavior is found in the critical line between the two polymerized 
phases: for x constant and slightly below the tricritical value, as x ls increased one passes 
from the Pb phase to the P phase and them back to the Pb phase. It may be understood 
that for high values of x the bilayer polymerized phase is favored, since in it the internal 
energy is lower, as compared to the regular polymerized phase, which is characterized by a 
larger entropy. 



V. CONCLUSION 



As shown above, the phase diagram of polymers placed on a bilayer with first neighbor 
interactions between bonds in different layers displays a variety of critical surfaces and 
multicritical lines and points, due to the fact that one non-polymerized and four distinct 
polymerized phases appear. Although to our knowledge such a system has not yet be studied 
experimentally, it is possible that an experimental situation similar to the model we study 
here may be found, with polymers confined to a region close to the interfaces between 
immiscible liquids, for example. 

As discussed in the introduction, magnetic bilayers have been much studied in the litera- 
ture by a variety of theoretical methods, including pairs of Bethe lattices similar to what is 
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FIG. 10: (Color online) Phase diagram for x\ = X2 = x and q = 3. Two tricritical points (orange) 
and a critical end point (blue) shown. 

done above [3] , and it is well known that polymer models may be mapped onto ferromagnetic 
n-vector models in the limit n — > (this limit is described in detail by Wheeler and Pfeuty 



17|). This correspondence may be generalized along the lines proposed for polymerization 



in poor solvents [18| , and in a particular limit may map into an effective n — > 0-vector model 
with higher order spin interactions between spins in distinct layers, and therefore it may be 
possible that similar phase diagrams to the ones found here may appear in related magnetic 
models, even for other values of n. Thus, there exists a possibility that phase diagrams 
similar to the ones found here may appear in appropriate magnetic bilayers. 

We notice that in all phase diagrams, with the exception of the first one shown in figure [3J, 
the most interesting features, such as multicritical loci, appear in the region of parameters 
where the interaction between bonds is attractive (x > 1), and this corresponds to the 
effective models for polymers in poor solutions mentioned in the introduction. 
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